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ABSTRACT 

00 

\ Ultra high energy protons accelerated at the shocks causing gamma ray bursts photoproduce 

^\ ' pions, and then neutrinos in situ. 1 consider here the sources of losses in this process, namely 

adiabatic and synchrotron losses by both pions and muons. When the shocks under consideration 
C I are external, i.e., those between the ejecta and the surrounding interstellar medium, I show that 

neutrinos produced by pion decay are unaffected by losses; those produced by muon decay, in 
the strongly beamed emission required by afterglow observations of GRB 971214, are limited in 
energy, but still exceed 10^^ eV. In particular, this means that ultra high energy neutrinos will 
be produced through afterglows. 
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1. Introduction 

The discovery of the afterglows of Gamma Ray Bursts (GRBs; Costa et ai, 1997; van Paradijs et ai, 
^ I 1997; Frail et al, 1997) and the disappearance of flares in their radio flux (Goodman 1997) have provided 

strong evidence in favor of the flreball model (Rees and Meszaros 1992), whereby mass is ejected from an 
as yet unknown source with large Lorenz factors, 7 ~ 100. Several authors have pointed out that in the 
extremely energetic events giving rise to these hyperrelativistic flows, high energy neutrinos are likely to be 
■ produced. A first likely source is the merger of the two neutron stars which may (Narayan, Paczynski and 

Piran 1992) give rise to the bursts themselves, as originally suggested by Eichler et al. (1989). Although 
insufficient to power the burst itself (Janka and Ruffert 1996), this may yet provide a copious source of 
low energy neutrinos, « 1 MeV. Another source of neutrinos comes from p — p collisions inside the shocks 
which give rise to the observed burst proper (Paczynski and Xu 1994); although it seems unlikely that 
these collisions may give rise to the burst itself, yet they may originate a flux of higher energy neutrinos, 
« 30 GeV, which is however currently unobservable (Ostrowski and Zdziarski 1995). 

A third, distinct source of high energy neutrinos exists. It has recently been pointed out that gamma 
ray bursts may be responsible for the acceleration of the highest energy cosmic rays observed so far (Vietri 
1995, Waxman 1995). It has been suggested (Waxman and Bahcall 1997) that these UHECRs may produce 
neutrinos through pion photoproduction, i.e., through the reaction 

p + 7 n + 

7r+ ^ + (1) 
^+ e+ + i^e + i^u , 



- 2 - 



inside the burst itself, i.e., against the 7-ray photons of the burst. While the probability of this happening 
is not large (Vietri 1995), yet Waxman and Bahcall (1997) pointed out that significant amounts of neutrinos 
with energy w 10^^ eV ought to be produced, and might be revealed by detectors of the AMANDA class. 
This detection is however made difficult by the fact that the Earth becomes opaque for this neutrino energy, 
so that the technique employing upward-moving muons cannot be used. But, through the same mechanism, 
and simultaneously with these lower energy brethrens, neutrinos in the energy range ~ 10^^ eV (UHENs = 
Ultra High Energy Neutrinos) ought to be produced; in Vietri (1998, paper I from now on) I showed that 
a significant fraction (/,r ~ 0.01) of all energy initially released in UHECRs ought to be lost to neutrinos 
emitted during the burst proper, or during the first week of the burst afterglow. The interesting point of 
this computation is that it appears that the next generation of satellite-borne detectors of large showers, 
such as AIRWATCH (Linsley 1997) may cover an area large enough to allow detection of UHENs within 
the first year or so of operation. 

Shortly before paper I was accepted, Rachen and Meszaros (1998) have presented a detailed analysis 
of losses limiting neutrinos' energies, in the internal shock scenario for GRBs, which applies for instance to 
the model of Waxman and Bahcall (1997). They conclude that these losses severely limit both the energies 
of individual neutrinos, and the total energy released by gamma ray bursts. Such losses arise because both 
the photoproduced pions, and the muons generated by the pion decay, before decaying suffer adiabatic 
and synchrotron losses, thusly imparting to neutrinos less than the energy they started out with. Given 
the short rest lifetimes of pions and muons, one maybe inclined to think these losses at all negligible, but, 
because of relativistic time-dilation, pions and muons moving with a Lorenz factor typical of ultra high 
energy cosmic rays (w 10^°) may survive more than 10* s, and cover sizeable distances (w 10^* cm) in the 
meantime. 

It is the purpose of this Letter to carry out the analysis of losses for external shocks, both during the 
burst proper and in the afterglow phase. It will be concluded that in the relatively tame environments 
generated by external shocks, neutrinos produced by pion decay are unaffected by losses, while neutrinos 
produced by muon decay are limited in energy by these processes, but still manage to exceed 10^® eV. The 
major departure from the discussion in paper I is due to the sensational discovery (Kulkarni et al., 1998) of 
the redshift z = 3.43 for the burst GRB 971214. This discovery forces us to choose different scaling values 
for bursts, for the following reason. For cosmological parameters f2 = and Hq = 50 km Mpc~^, the 
luminosity distance to this redshift is 6.2 x 10^* cm; for f7 = 1 it would be 8.4 x 10^* cm. The burst as 
observed by the BeppoSAX GRBM had a peak luminosity of 6 x 10"'' erg cm~'^, and a total duration 
of » 30 s (Heise et al., 1997), corresponding to a fluence of « 10~^ erg cm~'^. With the lower, above-stated 
luminosity distance this converts to a total energy release of 4 x 10^^ erg, roughly the binding energy of two 
merging neutron stars. While strictly speaking still not an impossibility, these values seem to imply some 
beaming, a natural consequence of many GRBs' models. In fact, assuming a beaming factor 9 f« 0.1, the 
total energy release becomes a more mundane 10^^ erg. Similar, although less extreme values have been 
derived, less cogently, for the burst GRB 970508 (Panaitescu, Meszaros and Rees 1998). These values will 
be adopted here. 



2. On the acceleration of protons in fireballs 

Before considering losses, I need to establish some properties of the environment, within the fireball 
theory of gamma ray bursts. The maximum proton energy, in the observer frame, was shown to be (Vietri 
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1995) 

,ma. = 2 X 10^° eV^l'^El'M^ezTe" (2) 

where the ejecta shell has, at the moment of impacting with the interstellar medium of density 
p = niTOp g cm~^, a Lorenz factor rj = 100?72, the shell explosion energy is E^i x 10^^ erg , 6 = O.lB^i 
the beaming angle, and ^ is the magnetic field energy density in units of the equipartition magnetic field 
energy density. As mentioned in paper I, this corrects a small error of Vietri (1995). When seen from the 
shell, this highest energy proton will have a Lorenz factor jp given by 

7p - 2 X 10%'^' Eli\Y'9-J/'e^' . (3) 

The two equations above correspond to Equations 29 and 28, respectively, of Vietri (1995). It is perhaps 
useful to remind the reader that they were obtained under the assumption that the major limiting factor of 
the protons' energies comes from the finite shell thickness: whenever the proton energy increases so much 
that g times its Larmor radius tl exceeds the shell thickness, the proton will cross the whole shell without 
being deflected backwards, and the acceleration cycle stops. Here, the magnetic field is assumed to be a 
fraction ^^/^ of its equipartition value, and the value g « 40 is taken from the numerical simulations of 
Quenby and Lieu (1989). 

However, at the referee's prompting, I can also show that the acceleration time-scale is shorter than the 
shell light-crossing time, which is the time-scale on which adiabatic losses set in, and also the time-scale 
on which the highest energy protons cross the shell before being deflected backwards. By establishing this, 
I shall in fact demonstrate that the two above equations correctly describe the highest proton energies 
achievable. The acceleration time-scale at relativistic shocks tr is shorter by a factor Q « 13.5 than the 
traditional acceleration time-scale at non-relativistic shocks (Quenby and Lieu 1989); thus 

Here matter on the two sides of the shocks moves with speeds Vi , V2 with respect to the stationary shock, 
A = grj^ is the typical deflection length, pi is the injection momentum, and pf the highest achievable 
momentum in the shock frame. Since A — gr^ — gcp/eB, I find 

Q{Vi^V2)Jp^ eB\Vi 1/2/ {Vi~V2)QeB\Vi ¥2)' 

Specializing now to a relativistic flow, for which Vi,V2 ~ c, and also Vi — V2 = qc ~ c, 1 obtain 

U^^. (6) 

On the other hand, for the highest energy protons, we have set, in deriving Eq. |^, A = gr^ = gcpf/eB ~ rtf^, 
with rth being the shell thickness; inserting this into the above equation, I find 

t. = ^«0.1^, (7) 
qCjc c 

which is as long as the shell light crossing time. So the acceleration is prompt enough to propel protons to 
the highest energies computed in Eq. H, within the time-scale rth/c, which is approximately the whole time 
the proton spends inside the shell, before reaching an energy so large that it cannot be turned backwards. 
Furthermore, proton adiabatic losses also become appreciable on the shell expansion timescale (Rachen and 
Meszaros 1998) rth/c, longer than the acceleration timescale Eq. 0. 
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It should be noticed that this result is independent of the specific values of g assumed; nor does it 
depend upon the magnetic field having achieved or not its equipartition value, ft is only the highest proton 
energy which depends on the parameters g and ^. Nor does it depend explicitly upon the shock Lorenz 
factor, provided the shock is relativistic of course. This is especially important when one considers the 
application of this same idea to the bursts' afterglows: the maximum proton energy is again provided by the 
finite shell thickness, which however, following Eq. ^, changes as the shell slows down. The reason for the 
continuing applicability of Eq. ^ is that, even during the afterglow, the shell thickness is given, in the shell 
frame, by rth — r/c where r is the current shock position, exactly as during the burst proper (Meszaros, 
Laguna and Rees 1994). This result can be phrased as follows: as the shock decelerates, the acceleration 
time-scale becomes longer, but so does also the shell crossing time and, provided the magnetic field builds 
up to the same, constant, fraction of the equipartition value ^, the two time-scales lengthen by the same 
amount, thus keeping the result of Eq. |^ still valid. 

The only hidden hypothesis in this argument is that appreciable magnetic fields exist on both sides of 
the shock to deflect the cosmic rays; in Vietri (1995) I suggested that this occurs when two or more shells 
collide. I believe the same suggestion is still valid for afterglows; I shall however discuss the implications of 
these multiple collisions in a different paper. 

Lastly I should remark that the results discussed in this section do depend upon the correctness of the 
assumption made by Quenby and Lieu (1989) that scattering of high energy protons at relativistic shocks 
is isotropic, and large-pitch angle. Since this provides a satisfactory explanation for the acceleration of 
cosmic rays at AGNs (Quenby and Lieu 1989), and it is borne out by theoretical arguments and simulations 
of interplanetary shocks (Moussas, Quenby and Valdes-Galicia 1982, 1987), it seems quite a plausible 
assumption to make. 



3. On fireballs 



The Equations ^ and ^ are however in an inconvenient form, the reason being that the parameter ni 
is unobservable. It is expedient to substitute for it the burst duration T — T2 x 100 s, a directly observed 
quantity which is related to it by the following argument. The shock with the interstellar medium forms at 
a distance Vgh from the explosion site, where (Meszaros, Laguna and Rees 1993) 

/ \ 1/3 

while the burst duration is of order 

T^^. (9) 

This is an old argument, originally due to Ruderman (1975), and more recently discussed by Sari and Piran 
(1997). Eliminating from the above I find 

and 

Tsh = irfcT = 6 X 10^'' cm T2 r^l . (11) 
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The magnetic field behind the shock, in the shell frame, is often assumed (Meszaros, Laguna and Rees 
1993) to be some fixed fraction S}^"^ of its equipartion value B^q, and is given by 

B = (87rCnimpc2)i/2^ = 54 G sH^i^^^'^T^'^'^ezlS}'^ , (12) 
where I have of course used Eq. ^ It is also convenient to rewrite Eq.|[ again using Eq. 10, as 



and Eq. ^ as 



7p = 3 X 10^ Eli'^^'T-'/'ezie^' 



3 X 10^0 eVEli'r^-'T-'/'ezie^^ 



(13) 
(14) 



4. Losses 



I have established above that protons suffer no adiabatic losses during their acceleration; in Vietri 
(1995) I had shown that also synchrotron and photopion losses are negligible. In Vietri (1998) 1 showed 
that the small photopion production leads, through the decays of Eq. 1, to the production of a small, but 
significant flux of prompt neutrinos. However, Rachen and Meszaros (1998) have shown that one ought to 
consider carefully whether synchrotron and adiabatic losses by muons and pions before their decay may 
significantly limit neutrinos' energies. 

The most significant losses will be due to muons, because of their longer lifetime (r^ = 2.2 x 10"^ s, 
= 2.6 X 10~* s). The maximum muon Lorenz factor 7^ before significant adiabatic losses set in is given 
by the following argument. A muon moving with Lorenz factor 7 in the shell frame appears to have a 
lifetime 7r^, because of relativistic time dilation. Within this time, the shell is expanding, thusly allowing 
the comoving magnetic field to decrease on a time-scale « rth/rth- The shell typical size rth in the shell 
frame is r/rj (Meszaros, Laguna and Rees 1993), where r is the shock instantaneous distance from the origin 



of the explosion; at the moment of the burst r = r^h (Eq. 11). The shell expansion velocity fth ~ c. The 



limiting Lorenz factor 7^ is found by equating these two time-scales: 

7m = ^ = 9 X 10"r2772 . (15) 

The condition that neutrino production is not affected, 7p < 7^, can be expressed as a requirement on the 
burst parameters: I obtain 

m > S.lEli^T-^^^ezl^'^^^^'^ . (16) 
For pion adiabatic losses, the only change comes in because of its shorter lifetime: 

7,2 > 0.5E}fT^'^^9Zy'e^' ■ (17) 

Synchrotron losses by muons are related to those of protons as follows. Calling 
ts — 1 yr (10"'^^/7)(1 G/BY the lifetime against synchrotron losses for a proton of Lorenz factor 7, 
the largest muon Lorenz factor 7^ before significant losses set in is given by (Rachen and Meszaros 1998) 

lis(^X = zlr,. (18) 



m 
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Here the muon mass is « OArUp. I find, using Eq. |l^, 



= 4 X IQi^i^ = 7.4 X 108i?5Y/'%'T|/'^_ir'/' • (19) 
B 



Again, the condition that these losses are irrelevant, jp < 7^, can be rewritten as a requirement on the 
burst parameters: I find 

ri2 > l.SEK^T^^^^ezl^^S,^^^ ■ (20) 

For pions the relevant criterion is 

772 > OAEli^'T-^^^ezTe^^ . (21) 

It is convenient to remark here that the use of 77 and T as independent burst parameters allows 
discussion of afterglow with the same formulae. During afterglows, the quantity T loses meaning because 
the afterglow emission is of course continuous. However, let us still define, by analogy with Eq. ^, T = 
where r and 77 are the instantaneous shock position and Lorenz factor. During the afterglow, for adiabatic 
expansion the shell Lorenz factor decreases as 77 cx r~^^^. Thus, the variation of T parametrizes the shell 
radius. I find that 

rp \ -3/8 



ri = m[j.J (22) 

describes the further evolution of the shell. Here 77^ and Ti are the initial values of the afterglow, at the 
moment of the formation of the external shock. This equation can also be established, for an adiabatic 
expansion in a constant density environment, by imposing the constancy through the afterglow of Eq. It 
is also identical to the time evolution law for adiabatic expansion with respect to observer's time, so that 
the fictitious quantity T can also be identified with Earth time, a physical quantity. In other words, by 
considering arbitrary values of 77 and T, I am considering points which do not really model any known burst, 
but which parametrize correctly later, afterglow moments of realistic bursts; the realistic initial models are 
linked to later values of rj and T by Eq. |2^. 

First I consider neutrinos produced directly by charged pion decay. In Fig. 1 I show the constraints Eq. 
|l7| and ^ for i?5i = 9_i = ^ = 1, plotted as solid lines. The dashed lines represent the afterglow evolution 
tracks, Eq. |2^. The figure covers roughly the first three afterglow hours. The dotted lines represent the loci 
of points of constant highest neutrino energies: when losses are negligible, these are computed multiplying 
highestproton energies, Eq. |l^, by 0.05. Neutrinos produced directly by charged pion decay will not be 
limited by the rather more stringent requirements on muon losses, Eq. This means that, should a burst 
fail to produce high energy neutrinos from muon decay, its neutrino flux will only be decreased (roughly) 
by a factor of 3 because the two neutrinos from the muon decay will be lost, at the same energy level. It 
can easily be seen from Fig. 1 that neutrinos of energy as large as 10^^ eV may be produced. 

Muon-produced neutrinos will be limited in energy both by synchrotron losses and by adiabatic losses, 
depending exactly on the model parameter values. The cross-over occurs when Eq. ^ equals Eq. 

mcr = O.SEli%-'/'9-J/'e^' (23) 
with adiabatic losses dominating whenever 772 > 772cr- 



When synchrotron losses dominate we see by comparing Eq. |19| with Eq. 22 that afterglow evolution 
occurs along lines of (roughly !) constant 7^ in the shell frame, so that the maximum achievable energy 
in the observer frame is given by 100 x 772 times Eq. Multiplying times the typical energy of outgoing 
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neutrinos, ~ 50 MeV we find that the neutrino energy is 

= 3.7 X lO^s eVE^^^/^rilT^^^e^iC'^^ ■ (24) 

Taking into account that during afterglow the shell Lorenz factor varies according to Eq. |2^, the typical 
neutrino energy is given by inserting this into the above equation: 

= 3.7 X 10^^ eVT^j^T^/^C'^^ (25) 

where of course r]2i and are the initial Lorenz factor and burst duration in units of 100 and 100 s, 
respectively, and the result is independent of post-burst time. Since the dependence on the initial parameter 
is so steep, it is enough that r]2i exceeds unity by a factor of 2 for to exceed 10^^ eV. When adiabatic 
losses dominate, the limiting muon Lorenz factor in the observer frame is given by 7^°''^ = 9 X 109 T2vl 
The locus of constant 7/1°'"^ is thus parallel to the top solid line of Fig. 1 (Eq. |l^), and to Eq. This 
gives a typical neutrino energy of 



e 



(""^l = 4.5 X 10" eVT2T]l . (26) 



However, afterglow evolution (Eq. |2^) is shallower than either of these; thus during the afterglow the energy 
of neutrinos produced becomes higher; this can be checked by computing the highest 7^°'"*'' for the afterglow 
evolution, Eq. I find 7/°''*'' oc T^/^. The normalization is such that the highest neutrino energy occurs 
on the top solid line of Fig. 1: inserting Eq. 16 into Eq. |lj, and considering that typical neutrinos carry 
away 0.05 of the proton energy I find 

= 5 X 10^8 eV EU'eZ^e^' ■ (27) 



In Fig. 2 I show as dashed lines tracks of constant energy for neutrinos produced by muon decay only, 
in the observer frame, together with the previously defined afterglow tracks (Eq. |2^, again dashed lines). 
The bottom solid line, Eq. |2^, marks the boundary between the region where adiabatic losses (above) and 
that where synchrotron losses dominate (below) . The top solid line marks the region (below the line) where 
adiabatic muon losses limit neutrinos' energies. At fixed T2, the neutrino energy is not monotonic with j]2 
because, above the region where adiabatic losses operate, each neutrino carries away a fraction « 0.05 of 
the emitting proton energy, which is a decreasing function of 772, for fixed ni (Eq. |l^). In either case, i.e., 
synchrotron or adiabatic losses, we see that it is still possible to produce high energy neutrinos, of order 
10^9 eV. 

In paper I, I showed that the fraction /^r of the photon luminosity Ly that goes into neutrinos is 

,300 keV 1 , , 

U = 0.01772-4 T2' (28) 

where « 300 — 10'^ keV is the typical burst observed spectral turnover energy. The important parameter 

and its time-dependence during the afterglow are not currently observed, nor are they in any way 
predicted by theory. A simple argument allowed to state that, if scales with the shell Lorenz factor 77 as 
oc 77" with a > 1, then remains constant or grows through the afterglow; if remains constant, then 
the neutrino luminosity — Jt^L^ oc t^^, the last proportionality being observed (Fruchter et al., 1997). 
Thus, even the bursts which, initially, suffer important muon losses, during the afterglow become effective 
neutrino emitters and, given that the neutrino luminosity scales only with t~^, the total neutrino fluence 
will not be much reduced by losses. The reason why I find now that, initially, losses are important is that, 
when I assume a significant amount of beaming, the maximum proton Lorenz factor, Eq. O, is greatly 
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increased, while the maximum muon Lorenz factor, Eq. [l^ for adiabatic losses, or Eq. |19| for synchrotron 
losses, is reduced or at most remains constant; it can be seen that in the isotropic case assumed in paper I 
the requirement on losses would have been less severe; however, as a partial compensation, for same total 
energy but significant beaming the maximum proton energy in the observer frame also increases, so that 
the end result is that, losses or not, neutrinos with energies exceeding 10^* eV are emitted. 



5. Discussion 

The propagation of ultra high energy cosmic rays in intergalactic space is significantly limited by energy 
degradation via pion photoproduction off the CMBR photons (Greisen 1966, Zatsepin and Kuz'min 1966). 
However, it has long been recognized (Wdowczyk, Tkaczyk and Wolfendale 1972) that observations of high 
energy neutrinos, thanks to their negligible cross sections for interaction with matter and photons, will 
eventually allow us to investigate the sources of cosmic rays way beyond this limit. With the development 
of Airwatch-class experiments (Linsley 1997), capable of monitoring from space atmospheric areas of order 
10^ krn^ in the search for extended air showers, and with detection efficiencies close to 1 for neutrinos (L. 
Scarsi, private communication), it seems the time for extending our search for sources of ultra high energy 
cosmic rays to the whole Universe has come. 

Two comments are in order. First, as already noted in paper I, the generic neutrinos (i.e., those 
produced by cosmic rays which managed to escape unscathed from their sources, but which produce 
pions off CMBR photons) should display dipole and quadrupole moments which, if UHENs come from 
GRBs, should be totally negligible (Fishman and Meegan 1995). This by itself may help rule out/establish 
an important alternative idea, that UHECRs have an origin in the Local Supercluster (Stanev et ai, 
1995). Second, the expected neutrino fluxes computed by Yoshida and Teshima (1993) for cosmologically 
distributed sources, leading to expected event rates for Airwatch-class experiments of 200 yr'^ , may be a 
significant underestimate. In fact, they assumed a distribution of sources of UHECRs arbitrarily limited to 
a redshift Zmax = 2. However, the recent, sensational identification of a redshift z = 3.43 for the burst GRB 
971214 (Kulkarni et ai, 1998) makes it clear that significant fluxes of ultra high energy neutrinos may be 
at very large redshifts, thus not violating the constraints on fluxes of UHECRs observed at Earth. Thus 
the building of Airwatch-type experiments with significant neutrino-detection capabilities becomes an even 
more exciting prospect. 

In this paper, I have strengthened the point made in paper I that a fraction /jr ~ 0.01 of all ultra high 
energy neutrinos should correlate (within about a week) with the burst where it was emitted, by showing in 
some detail that losses do not inhibit the production of high energy {'^ 10^^ eV) neutrinos, simultaneously 
with the burst or its afterglow. A search for ultra high energy neutrinos ought to be a signiflcant test of the 
hypothesis that UHECRs are generated in GRBs. 

It is a pleasure to acknowledge helpful comments from J. Rachen and P. Meszaros, and from an 
anonymous referee. 
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Fig. 1. — Typical energies for neutrinos produced by charged pion decay, as a function of bursts' parameters. 
Top solid line (Eq. ^7|): pion adiabatic losses are negligible in the upper region. Bottom solid line: Eq. 
pT| . Below this line synchrotron losses become important. Dashed lines: representative tracks for afterglow 
evolution, Eq. Dotted lines: loci of points of constant neutrino energies in the observer's frame. The top 
dotted line corresponds to €1, ~ 1.5 x 10"'^* eV , while the bottom dotted line to tu = 1.5 x 10^^ eV . These 
values are obtained by multiplying Eq. |lj times 0.05, the typical yield for the neutrino energy. The figure 
covers, roughly, the first three afterglow hours. 
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Fig. 2. — Typical energies for neutrinos produced by muon decay as a function of bursts'parameters. The 
bottom solid line (Eq. p3|) separates the region where adiabatic losses prevail (above the line) over synchrotron 
losses. The top solid line solid line (Eq. |l^ separates regions where adiabatic losses are important (below 
the line) from that where where they are negligible (above the line). Dotted lines: loci of constant neutrino 
energy, in the observer's frame. The top dotted line corresponds to = 1.5 x 10^^ eV, while the bottom 
dotted line to e,y = 4.5 x 10^^ eV. Dashed lines: representative afterglow tracks, Eq. Since afterglow 
tracks are shallower than constant neutrino energy tracks, it can be seen that later in the afterglow emitted 
neutrinos will have larger energies. 



